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AND A FAMILY OF ALGEBRAIC CURVES 
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Abstract. We present a set of conditions which, if satisfied, provide for a complete asymptotic analysis of 
random matrices with source term containing two distinct eigenvalues. These conditions arc shown to be 
equivalent to the existence of a particular algebraic curve. For the case of a quartic external field, the curve 
in question is proven to exist, yielding precise asymptotic information about the limiting mean density of 
eigenvalues, as well as bulk and edge universality. 



1. Introduction 

1.1. Random Matrices with Source. Consider the probability measure 

(1) fx n (dM) = —e~ n Tt(v{M)-AM) dM; 

defined on n x n Hermitian matrices, M. Here dM denotes Lebesgue measure on the matrix entries, 
dM = rij=i dMjj Oi<j<fe<n dHe(Mjk)dhn(Mjk) > the matrix A is a fixed n x n matrix, and the parameter 
Z n is a normalization constant chosen so that Q is a probability measure. 

This defines an ensemble of random matrices, in which the matrix A plays the role of an external source. 
The function V should be real and grow sufficiently rapidly that the above measure possesses all finite 
moments; it is typically assumed to be a polynomial. 

The family of measures described by Q are not invariant under unitary transformations, and the anal- 
ysis of eigenvalue statistics under these measures is in its infancy in comparison to the so-called "unitary 
ensembles" for which a great deal is known. 

Following the work of Brezin-Hikami , [Hi , and P. Zinn- Justin > a team of researchers ( Bleher 

and Kuijlaars [HJ,01 and Aptckarcv, Bleher and Kuijlaars PP) have considered the large n behavior of 
eigenvalue statistics under from the point of view of Riemann-Hilbert problems. 

In 1 1 3j , Zinn- Justin showed that the eigenvalues of such random matrices are in fact a dctcrminantal point 
process, and in |2j, the authors showed that this representation may be described in terms of the solution of 
a matrix Riemann-Hilbert problem, the size of which depends on the number of distinct eigenvalues of the 
matrix A. 

In the present paper we shall assume that the matrix A possesses two eigenvalues, —a and a, with 
multiplicities ri\ and ri2, respectively (with n\ + = n). In this case, the Riemann-Hilbert problem is as 
follows: 

Riemann-Hilbert Problem 1.1. (a) A is analytic on C \ R. 
(b) The boundary values of A satisfy 

(2) A+{x) = A_(x 
where 

nVi p — n(V(x)+ax) 



(3) W! = e~ nVl = e 
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(c) As z — ► oo, we have 

(4) ^ Z ' = ( / + °(J)) |° ^ ° j' 

In terms of this Riemann-Hilbert problem, the probability measure on eigenvalues induced by (y| may 
be re-written as follows: 



(5) dfi n (x!, . . . ,x n ) = det (Knixi^Xj))^^^ d n x, 

[T5| . [2] where now 

-n{V(x) + V(y)) f 

*»(*■») = 2m{x _ y) {e^[Y(y)-iY( X )} 21 



+e" 



**» [y( y )- 1 y(x)] 31 } 



The formula Q is not just a concise representation of the probability measure, it turns out that all 
statistical properties of eigenvalues can be related to K n . For example, 

Prob {no eigenvalues in (a, b)} = det (1 — ^n)i,2u a mi 

where IC n is the integral operator with kernel K n (x, y): 

fcnf = / K n (x,y)f(y)dy. 

J a 

The representation of the correlation functions in terms of the kernel K n is due to Zinn- Justin , and the 
representation of this kernel in terms of a Riemann-Hilbert problem is due to Blchcr and Kuijlaars 2 . 

Remark: In the case that the matrix A possesses p distinct eigenvalues with multiplicities m,ri2, • • • , n p , 
the above representations generalize directly, but the associated Riemann-Hilbert problem isp+1 x p + 1. 

I n E3- the authors considered the case that V is a quadratic. In the series of papers 0, and jlj, 
the authors also considered the Gaussian case, but from the point of view of Riemann-Hilbert problems. 
The goal of these works was to study the behavior of the eigenvalue statistics in the large n limit. In |14j . 
Zinn-Justin also studied the large n limit, for the case of general V. In that work, one important issue 
was the description of the limiting density of states in terms of a function for which an existence theorem 
is lacking. He explains very carefully the analyticity properties and branch-cut structure required of this 
function, in the so-called "one-cut" case, under the assumption that the density of the eigenvalues of A has 
a smooth limit when n tends to oo. 

Here we consider the case that A possesses two eigenvalues, of multiplicities n\ and ri2, and our interest 
is in the behavior when n, m, n,2 — * oo so that — — > xj. The main goal of our work is to provide a new and 
explicit characterization of the limiting density of states in terms of an algebraic curve. The following is a 
summary of the results in this paper: 

(1) In Section|3 we present a list of conditions (referred to as an "ideal situation") which, if true, yield a 
transformation of the Riemann-Hilbert problem to a "normal form" , that is, a Riemann-Hilbert 
problem from which subsequent asymptotic analysis (for n — > oo) is straightforward. 

(2) In Sectional we show that if these conditions are satisfied, then there exists an algebraic curve whose 
roots yield the desired transformation. The curve is always of the form 

(6) w 3 -V'{z)w 2 +C 1 (z)w -Co(z) = 0, 

where C\ and C 2 arc analytic functions of z. 

(3) For arbitrary polynomial external fields, the algebraic curve may be determined up to a finite number 
of free parameters (see Section @J. 

(4) In Sectional we consider the Gaussian and Quartic cases. For the Gaussian case, this algebraic curve 
has no free parameters, and in the case n\ = n-i = n/2, is equivalent to the algebraic curve used by 
Blchcr and Kuijlaars 0E] and by Aptekarcv, Bleher and Kuijlaars [fl. 



(5) For the quartic case, V = x 4 /4, with ri\ = n-i = n/2, we prove that for a sufficiently large there 
exists a choice of the free parameters so that the algebraic curve yields the desired transformation. 
Each of the two relevant roots possesses a single branch cut, which is the "one-cut" case. 

(6) In Section \7\ in the "one-cut" case, we show that if the desired transformation exists bringing the 
Riemann-Hilbert problem into a normal form for subsequent Ricmann-Hilbert analysis, then the 
analysis of Blcher and Kuijlaars [3] can be applied directly. 

(7) It follows that if the desired transformation exists, then (i) the limiting mean density of eigenvalues 
exists and is described in terms of the algebraic curve, (ii) bulk universality holds true, and (hi) edge 
universality hold true (see Section |SJ. 

(8) The combined results of Sections |H1 and [7| show that for the quartic case V(x) = x A /A, with n even 
and n-i = n% = n/2, bulk and edge universality still hold true, for all a sufficiently large. 

Remark: In the analysis of singular limits of 2 x 2 Riemann-Hilbert problems associated to integrable 
systems and random matrix theory, there has often arisen an existence problem; a transformation is re- 
quired satisfying a number of nonlocal properties, in order to transform to a new Riemann-Hilbert problem 
from which subsequent asymptotic analysis is straightforward (see, for example, [J]). The existence of this 
transformation has been shown to be equivalent to the existence of an equation of the form 

F 2 + W(z)F+U(z) = 0, 

in which (i) the functions W and U are analytic and defined implicitly in terms of the desired transformation, 
(ii) all branch points lie on the contour of the Riemann-Hilbert problem, and (iii) the roots F±(z) of this 
equation satisfy a variety of inequalities and relations on various subsets of the contour. 

The proofs of existence of these transformations (or equations) are nontrivial, and require techniques 
ranging from WKB analysis of odes, through the theory of level trajectories of quadratic differentials, to the 
theory of logarithmic potentials and equilibrium measures. 

Viewed in this light, the present paper explains that for the analysis of 3 x 3 Ricmann-Hilbert problems, 
the existence of a suitable transformation is equivalent to the existence of an equation of the form JSJ) whose 
roots also satisfy a host of properties. For quartic potentials, we establish an existence theorem via algebraic 
geometry. 

2. LARGE n ASYMPTOTIC ANALYSIS OF THE RlEMANN-HlLBERT PROBLEM 11.11 

This paper studies the large n asymptotic behavior of the solution to the Ricmann Hubert problcm ll.il 
under the assumption that the functions Vj are entire. Furthermore, we will assume that the integers nj 
grow to oo with n, in such a way that 

(7) —^Xj, 3 = 1,2. 

n 

2.1. First step in the asymptotic analysis. In analogy with the more standard case of orthogonal 
polynomials (see, for example, 0-0), we seek a vector of functions 

(8) G(z) = (g (z), gi (z),g 2 (z)) 
satisfying the following set of assumptions. 

(1) There are 2 disjoint sets, I\ and I2, each a compact subset of R,with Sj denoting the supremum of 
each set. 

(2) For each j = 1, 2, the function gj is analytic in C \ (—00, Sj). 

(3) For j = 1,2, each function gj behaves as follows at 00: 

(9) g j (z)=logz + 0(z- 1 ). 

(4) The functions gj satisfy 

(10) 5o = . 9i +.92- 

In addition to these basic conditions, there will be a collection of important conditions on these functions, 
collectively referred to as the "Ideal Situation" . The goal is to deduce that there exists a unique vector G 
that can be used to transform the Ricmann-Hilbert problem 11.11 to a new Ricmann-Hilbert problem in a 
form suitable for asymptotic analysis. To that end, the transformation that we seek is as follows: 
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(11) B ~ diag ( e -^o/2 )e nili/2 )e n 2 i 2 /2^ A(z ) diag (^-n(g -l /2) > ^ ( 9l -h/2) _ g n 2 ( S2 _< a /2) 

where £ = (£oj — €1, —^2) is a vector of constants, which will be chosen along with the vector G. For now, 
the only restriction on these constants is that ti£q = n\l\ + n.2^2- 

The matrix B inherits jump relationships across R, that are determined by the Riemann-Hilbcrt problem 
11.11 for A, along with the boundary value behavior of the vector G. Therefore, B satisfies the following 
Riemann-Hilbcrt problem. 

Riemann-Hilbert Problem 2.1. (a) B is analytic on C \ R. 
(b) The boundary values of B satisfy 

/g— n (9o~ Bo) e n 9o +™i9i~ —nVi—nxii e ng„ +n 2 g} -riVi-n^l^ 
e ni(9?-Bl) 

1 e m{gt-si) 



(12) B+{x) = B^{x) 



(c) As z — > 00, we have 
(13) B(z)=[I + 



1 

z 

Note that for ease of notation, we have introduced the constants tj, j = 1,2 (without hats). These are 
defined by l x = «%fcm4, and l 2 = " io +" 2 ^ 2 . 

J 1 2ni 1 ^ 2ri2 

Ideal Situation 

In analogy with the case of orthogonal polynomials, a detailed asymptotic analysis would be possible if 
the following were true. 

A. (As mentioned before) There are 2 disjoint sets, I\ and I2, each a compact subset of R, so that for 
each j = 1, 2, gj is analytic on C \ (—00, 5j) (where Sj is the supremum of the set Ij). 

B. The function gi and its derivative g[ are continuous in the closure of C\ (—00, 5%), and the functions 
<72 and g' 2 are continuous in the closure of C \ (—00, 62)- That is, gj and g'j achieve their boundary 
values in the sense of continuous functions. 

C. For j — 1, 2, we have that for z £ Ij, the following two properties hold true: (i) 4 [g^ — gj) is real, 
and decreasing as one traverses Ij according to its orientation, and (ii) 

(14) ngo + n g+ - nVj - rijlj = 0. 

D. For j — 1, 2, we have that for z £ R \ Ij, the following two properties hold true: (i) i [g^ — gj) is 
constant, and (ii) 

(15) Re [ng^ + njgj — nVj — rijlj] < 0. 

Remark: The reader may wonder why conditions A-D above have been described as an ideal situation. 
The explanation requires considerable calculation, but the basic idea is this: on each subinterval Ij of R, the 
jump matrix may be factored in a manner quite analogous to the matrix factorization that has been used 
for asymptotic analysis of 2 x 2 Riemann-Hilbert problems in approximation theory over the past 8 years. 
The theory for the case of 2 x 2 matrices can be essentially embedded in the more general setting. The 
real evidence that this works is in the papers [Hj, [Hjj Q, and 0j, where the authors prove in special cases 
that the conditions A-D can be achieved, and then proceed to use this information to compute a complete 
asymptotic description for the Riemann-Hilbert problem. 

In Section we will carry out this Riemann-Hilbert analysis, for the "one-cut" case. Indeed, the reader 
will see that in this case, the subsequent Riemann-Hilbert analysis is virtually identical to the analysis 
carried out by Bleher and Kuijlaars [Hj- 

Remark: The "ideal situation" described in A-D above does not always hold true. Even in the Gaussian 
case, for a < 1, it turns out that the desired sets I\ and I2 are not subsets of R, but rather they are pieccwise 
analytic arcs (actually each is a union of two line segments). The more general setting can be described as 
follows 



• The Riemann-Hilbert problem 1 1.1 1 is modified so that the jump is not across the real axis, but rather 
is across two contours, Ti and IV 

• The ideal situation would be if the sets Ij , no longer subintervals of R, become unions of picccwisc 
analytic arcs, subsets of the contours Ti and T 2 , so that the conditions described in A-D above are 
true, on the sets Ij. (The conditions in B must be generalized so that the functions gj are analytic 
with continuously diffcrentiablc boundary values off semi-infinite subsets of the contours T\ and T2-) 

For V convex and a sufficiently large, it seems to be the case that the more down-to-earth conditions 
of A-D, with most action taking place on the real axis, should capture the situation, but this remains to 
be seen. On the other hand, when the more general setting does occur, the existence of an algebraic curve 
may be deduced as well, and the actual form of the equation depends only on the external field V, and the 
parameters a, n\ and n-i. 



Using the fact that ngo = n\g\ + 71-232, we may rewrite the condition (|14|) as 

(16) rij [gj + gf) + n k g k - nVj - n,j£j =0, z € Ij, kj = 1,2, k^ j. 

In the next Section we will show that the collection of boundary relations contained in (|16|) are intimately 
related to an equation of the form 

(17) w 3 + ci2(z)w 2 + a\{z)w + ao(z) = 

in which the coefficients cij(z) are analytic for all zsC. We will show that the roots of this polynomial are 
expressible explicitly in terms of the c^-'s and V^'s. 

Remark: For the case of p x p Riemann-Hilbert problems associated to Random Matrices with Source, 
the collection of boundary relations analogous to l|16|) should be related to an equation of the form 

(18) w p + a p - 1 (z)w p - 1 + ■ ■ ■ + a (z) = 0. 

3. An algebraic curve 
We begin by defining /1 and f 2 as follows: 

77 ■ 

(19) fj(z) = ^g'j(z), j = 1,2. 

n J 

The boundary relations (|16f) yield, upon differentiation, that f\ and /2 satisfy 

(20) (/f + /+) + h - V{ = 0, z€l u 

(21) (/ 2 " + /2 + )+/i-^' = o, zel 3 . 

In addition, if the condition that \ (g^ 1- — gj) is constant on R \ Ij holds true (see part D of the "ideal 
situation" of the previous section) then clearly for each j = 1, 2, fj is analytic in R\ Ij. 

Lemma 3.1. Assuming the hypotheses A-D set forth in the "ideal situation" of the previous section, we 
conclude that the functions fj satisfy the following relations for all z £ C 

(22) h(zf + h{z)h{z) + h{zf - h{z)V{{z) - h{z)V 2 '{z) = A(z) 

(23) h{z)h(z) (h(z) + f 2 {z)) - V{{z)h{z) (V{{z) - h(z)) - V 2 '(z)f 2 (z) (V 2 '(z) - h{z)) = B(z), 
where A(z) and B(z) are entire functions, implicitly defined as follows: 

(24) A( Z ) = J- / iim^l±M^M ds 

2ni J 1 s — z 

(25) B{z) = j_ r YMMil - + (jj(f) - h( s )) dc 

2iri J 7 s - z 

and in both \2J$ and \25)) . the contour of integration 7 may be taken to consist of a finite union of simple 
closed Jordan curves which do not intersect and whose interiors are disjoint, each oriented in the clockwise 
direction, one subset of these contours encircling the intervals comprising I\, a second subset encircling the 
intervals comprising I 2 , and one final contour encircling the point z. (Of course, since V\ and V2 are entire, 
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this may be deformed into one contour, a large circle oriented in the clockwise direction, which encircles the 
sets I\ and I 2 , and also encircles the point z.) 

Proof 

Let H 2 = fi(z) 2 + f 1 (z)f 2 (z) + f 2 (z) 2 -f 1 (z)V{(z)-f 2 (z)V^z). As defined, H 2 is analytic in C\(h U J 2 ). 
However, the jump of H 2 across either of these sets may be computed. Indeed, across I\, we have 

(26) H+(z)-H 2 -(z) = {f+(z)-f 1 -(z)){f+(z) + f 1 -(z) + Mz)-V;(z)} = 0. 

where in the last equality, we have used the boundary relation l|2U|) . In exactly the same way one may verify 
that H 2 also possesses no jump across I 2 . Since we know that f\ and f 2 are bounded, H 2 can have no 
isolated singularities, and so Morera's theorem tells us that H 2 is in fact entire. Equation IL'l'l) with A(z) 
defined in l|24|) follows by Cauchy's theorem. 

The proof of (|23|l goes along the same lines. One begins by defining 

(27) H 3 = h(z)Mz) (A(z) + f 2 (z)) - V{(z)A(z) (V{(z) - h{z)) - V 2 \z)f 2 {z) {V 2 \z) - f 2 {z)) , 

and observing that H3 also possesses no jumps across I\ U I 2 , and hence it must be entire. Then (|23|) with 
B(z) defined in (|25|) follows by Cauchy's theorem. □ 

Remark: The reader will note that the assumption that V\ and V 2 should be entire is for comfort; for 
example, if V{ and V 2 are meromorphic functions, then the only modification to the theorem is that A{z) 
and B{z) are also meromorphic, as they might have poles at those points that are poles of V{ and V 2 . The 
integral representations l|24[) and (|25|l for A and B remain true, with the contour of integration taken to 
consist of 3 simple closed Jordan curves that do not intersect each other, each oriented in the clockwise 
direction, one encircling I\, the second encircling I 2 , and the third encircling the point z. (It may not be 
possible to deform these into a single contour without introducing additional residues.) 

Remark: In order to conclude that H 2 can have no isolated singularities on I\ or I 2 , we don't actually 
need that the boundary values be bounded. The Lemma remains true if the boundary values of H 2 lie in 
L p (Ij) for some 1 < p < 2. This may be translated into conditions on the functions fj, but this is, at the 
present time, not particularly useful. 

Remark: We note in passing that not all of the assumptions of A-D are required for the proof of this 
Lemma, and in fact what is actually required is A, B, (|14f) . and item (i) of D. 

The following uniqueness result follows immediately from Lemma 13.11 

Theorem 3.1. Suppose that there are two disjoint sets, I\ and I 2 , each one bounded and consisting of a 
disjoint union of intervals. Then the boundary value problem for $20\) - $21\) (with, for j = 1,2, fj analytic in 
C\Ij, and the boundary values of fj being bounded) possesses at most one solution. 

Proof If there is another pair of functions f\ and f 2 that satisfy the same boundary value relations and 
asymptotics, then their difference Fj := fj — fj, j = 1, 2, satisfies the following relations. 



(28) F++Ff+F 2 = 0, zeh, 

(29) F++F-+F 1= 0, zel 2 , 

(30) Fj = ^,z^ 00. 
Thus, Lemma 13. II implies that for all z£C, 

(31) Fi + F1F2 + F| = 0, 

(32) F 1 F 2 {F 1 +F 2 )=0. 



These identities finally imply that Fj = for j = 1,2, and we have proven uniqueness. □ 

Lemma 13.11 clearly demonstrates that the boundary value relations H2Q(I - (|21(I are related to algebraic 
equations; it is useful to recast these relations in terms of the following result . 







Theorem 3.2. Define r\(z), r%(z), and r$(z) as follows: 




Then, under the assumptions of Lemma \S.l\ the following holds true: 

(36) E(w, z) := (w - ri) (w - r 2 ) {w — r 3 ) = w 3 - C2{z)w 2 + Ci(z)w — C (z), 



where 



(39) 



(37) 



(38) 



C 2 (z) = -(V((z) + Vi(z)) 

C 1 (z) = ^-(4A(z) + (V{(z)-Vi(z)) 2 ) 

C Q = l Ua(z) (V{(z) + Vi{z)) - 8B(z) - (V{(z) - Vi(z)f (V{(z) + V^(z)) 



Proof 



The proof is a straightforward exercise in algebra: one must compute the coefficients of each power of w 
in the expression for E, using the definition (|33|) ~ H35|) of rj, along with 12211 and l|23(l . So, for example, (|37|l 
may be verified as follows: 



from Q35[l . The relations (|38|l and 139|) may be verified in a similar fashion. □ 

Remark: As with Lemma f3. II the assumption that V\ and V2 should be entire is for convenience, and 
in particular if V\ and V% are meromorphic, then the Theorem remains true (but with A and B adjusted to 
admit the possibility of poles at the pole locations of V{ and V 2 '). 

Remark: It is important to observe that Theorem 13.21 yields nontrivial information about the functions 
<7i and 32 • Indeed, the following result provides a drastic reduction in the possible complexity of the sets I\ 
and 12- 

Theorem 3.3. Assuming the hypotheses A - D set forth in the "ideal situation" of the previous section, each 
of the sets Ii and I 2 consists of at most a finite union of intervals. 

Proof Since fi is analytic in a neighborhood of 00, and since I\ is a contour with finite arclength, f\ may 
be expressed as a Cauchy integral of its jumps across I\. The same holds true for f 2 , and so we have 



On the other hand, the functions r%, r 2 , and r 3 are the roots of a function E(z, w) that is a cubic polynomial 
in w with coefficients that are entire functions of z, and so its roots are analytic functions of z, with at worst 
a finite number of branching points in any finite part of C. From the definition of the r^'s, we learn that the 
functions f\ and f 2 can have at most a finite number of branching points in the entire plane. Each endpoint 
of a sub-arc of Ij clearly corresponds to a branching point, and so we conclude that I\ and I2 must each 
consist of at most a finite union of intervals. □ 

The conditions A through D above, when satisfied, should in principle allow one to carry out a complete 
asymptotic analysis of the associated multiple orthogonal polynomials, and obtain complete control of the 
asymptotic statistics of the eigenvalues. What has not been pinned down completely under these conditions 
is the behavior near the endpoints of I\ and I 2 ■ This edge behavior depends on the order of vanishing of the 
quantities ng$ + njg^ — nVj — nj£j and i (cj^ — gj^j at the endpoints. The following additional condition, 
expected to be the generic case, implies that all edge behaviors are governed by the Airy equation. We will 
take (a, (3) to be one of the intervals comprising the set (which is now taken to consist of a finite union of 
disjoint intervals). 



(40) 



C 2 (z) = r 1 +r 2 +r 3 = \(V{(z) + Vi(z)), 



(41) 
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E 1. The behavior near each left endpoint a is as follows. There is c a > so that 

• g${z) + Xjg~j{z) — Vj(z) — xjij = — c a \z — a| 3 / 2 + O (|z — a| 5 / 2 ) for z near a and z < a. 

• i (gj~(z) — g~(z)) = c a \z — a\ 3 ^ 2 + O (|z — a| 5 / 2 ) for z near a and z > a. 
E 2. The behavior near each right endpoint (3 is as follows. There is cp > so that 

• Qo(z) + xjg+(z) - Vj(z) - xjij = -c$\z - /3| 3/2 + (|z - /?| 5/2 ) for z near (3 and z > a. 

• i ( 5 +(z) - gj(z)) = c p \z - /?| 3 / 2 + O (\z - /3| 5 / 2 ) for z near /3 and z < (3. 



4. Polynomial external fields 



According to the notation and definitions of Sections 2 and 3, we have V{ = V'(x) + a and V 2 ' = V'(x) — a, 
and hence one may verify that the desired algebraic curve has the form 



(42) 



B] =0. 



w 3 - (V'(z)) w 2 - [A + a 2 ] w-[(A- a 2 ) V'(z) 

The quantities A and B are entire functions, as yet undetermined (aside from the fact that they satisfy (|24l) 
and EH ). 

Now if it is the case that V(x) is a polynomial, more information can be obtained about A and B. In this 
case, they satisfy 



(43) 
(44) 



A(z)= Poly [-V'(z) (fi(z)+f 2 (z))], 

B(z) = Poly f- (V'(z)f (h{z) + f 2 {z)) + V'(z) (/ 2 + 2a (f 2 - h) + f 2 ) 



where Poly(Q) refers to the polynomial part of the Laurent expansion of Q valid for z near oo. Clearly 
then, A is a polynomial of degree dcg(T^) — 2, and B is a polynomial of degree 2deg(T^) — 3. So, if V is a 
polynomial, then these considerations determine the algebraic curve in question, up to the choice of at most 
3deg(y) — 5 undetermined coefficients. 

Amazingly, the algebraic curve 1)42(1 is actually of much lower degree than might be expected. Indeed, 
some algebra shows that the curve has the simplified form 

(45) w 3 - (V\z))w 2 +d{z) w-C = 0, 

where C\(z) = Poly [V'(z) (/i(z) + f2(z))] — a 2 is a polynomial of degree deg(V^) — 2, and Co is a polynomial 
of degree deg(V) — 1. To see that this latter claim is true, note that one has the representation 



(46) 
(47) 



C (z) = -a 2 V\z) - Poly [V(z) (f 2 + 2a (/ 2 - h) + /f )] + 



Poly 



(V'(z) f (A(z) + /,(*)) - V'(z)Pdiy [V'(z) (A(z) + Hz))} , 



which is easily seen to be a polynomial of degree deg(V) — 1. 



5. The existence problem. 

We have shown that if there exists sets I\ and I2 and functions /1 and f 2 satisfying the conditions of 
Lemma |3.1I then the functions /1 and f 2 are related to the roots of a particular algebraic curve. On the 
other hand, for the purposes of asymptotic analysis of Ricmann-Hilbert problems, what is really required is 
the existence itself: given external fields V\ and V 2 , one needs to know that functions /1 and f 2 , together 
with sets I\ and I 2 , exist so that the original Ricmann-Hilbcrt problem 11.11 may be "regularized" into the 
form of Riemann Hilbcrt problem 12. II and subsequent asymptotic analysis may be carried out. 

If one knew the algebraic curve explicitly, then by careful analysis one could prove that (i) sets I\ and 
I 2 exist (they are the branch cuts of the algebraic functions determined by the algebraic equation), (ii) 
functions /1 and f 2 exist (they are the actual algebraic functions) and (iii) all the required reality, and 
positivity/negativity conditions are satisfied (this is actually the most difficult part of the analysis). 

Unfortunately, except in special Gaussian cases (see the next Section), the functions C\ and C 2 are not 
completely determined by the considerations leading to, nor by the proof of, Theorem 13.21 For example, if 
V\ and V 2 are polynomials, then one can prove (again, see the subsequent Section for concrete examples) 
that C\ and C 2 are polynomials in z, whose leading coefficients are explicitly known. However, the explicit 
determination of the lower order terms does not follow from these considerations. 
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In other words, one has not shown existence of sets Ij and functions fj unless one can prove that within 
the family of curves parametrized by the coefficients left undetermined by Theorem 13.21 there exists an 
algebraic curve which satisfies all the necessary conditions. We will formulate this as an existence problem. 

Existence Problem 5.1. Given polynomial (or more general) external fields V\ and V 2 , prove that there 
is an algebraic curve of the form $36)) whose roots yield functions f\ and f 2 as in $33)) - \3~lj\) so that items 
A-D of the "ideal situation" described near the end of Section 2 (or their generalization to the case where 
the curves in question are not necessarily real) are satisfied. 

6. Existence theorems for quadratic and quartic external fields 



The Gaussian case. 

The Gaussian case corresponds to V = x 2 /2. In this case, one may verify that 



(48) A(z) = -I, 

(49) B(z) = -z. 

And so, we see that the desired algebraic equation is 

(50) w 3 - zw 2 - (a 2 - l) w + a 2 z + a (2x 2 - 1) = 0. 



We observe that in the Gaussian case, there are no undetermined coefficients at all, and, it is straightforward 
to verify (see, for example, 00, for the case x\ — X2 = 1/2) that there is A so that for a > A, i|50|) possesses 
roots that yield functions gj and sets Ij so that the conditions A-D hold true. (In the case x\ = x 2 = 1/2, 
it is also known that for a < A, the equation (|5()[) determines two piecewise analytic arcs as described in the 
third remark following the description of the conditions A-D in Section [SHI])- 

The quartic case, with x\ = x 2 = 1/2. The quartic case corresponds to V(x) = x 4 /4. In this case, one 
may verify that 

(51) Ci = z 2 + C\Z + c , 

(52) C = -a 2 z 3 + c 2 z 2 + c 1 z + c . 

If we consider the special symmetric case in which x\ = X2 = 1/2, then there is an additional symmetry 
satisfied by the functions f\ and f 2 , namely that 

fi(z) = -M-z). 

(This can be seen by considering Fi(z) = — f 2 {— z) and F 2 (z) = —fi(—z), showing that F\ and F 2 satisfy 
the same boundary value relations H20() - (|21|l as /1 and f 2 , and then using Theorem 13. II ) 

This symmetry in turn implies that Co is an odd polynomial, and C\ is an even polynomial, and hence 
the algebraic curve becomes 

(53) w 3 - z 3 w 2 + (z 2 + a)w + a 2 z 3 + f3z = 0, 

and there are only two undetermined coefficients, a and (3. 

In order to arrive at a problem of the form of Riemann-Hilbert problem 12.11 the existence of sets I\ and 
I 2 , along with functions f\ and f 2 must be established, as described in Section 3. As a demonstration that 
this is possible, we have the following result. 

Theorem 6.1. There is A > so that for all a > A, there exists coefficients a(a) and /3(a) so that the 
following properties, pertaining to equation if 53)) with a = a(a) and (3 = (3(a), are true. 

(1) The three solutions {fj{z)} 3 =1 of the equation )53]) . determined uniquely as algebraic functions of 
z for z G C that possess standard meromorphic Laurent expansions at 00, satisfy the following 
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asymptotic behavior for z — > oo : 

(54) r 3 = z 3 -- + 0(z- 2 ) , 

z 

(55) r 2 = a+^+0(z~ 2 ), 

(56) n = -a + L + o (z~ 2 ) 

(2) The functions f\ and f 2 , defined by 

(57) h(z) = ri (z) + a, f 2 (z) = r 2 {z) 

possess the following asymptotic expansions for z — ► oo : 

1 

2~ 

(3) There are two real positive real numbers 71 = 71(a) and 72 = 72(a) so that the function f 2 may be 
taken to be analytic in C \ I 2 where I 2 = [71 , 72] . 

(4) The function f\ may be taken to be analytic in <C\I\ where I\ = [—72, — 71]. 

(5) The following boundary relations hold true: 

(59) (/r + /+) + f 2 - V( = 0, zel u 

(60) (f 2 + /+) + h - V 2 ' = 0, zel 2 . 

(6) For z G Ij, we have 

(61) i(/+(z)-/r( z ))<0. 



(58) f j = - +0 (z- 2 ), j = 1,2. 



(7) For z£l\Ji, we have 



(62) / i?e (/f ( S ) + /+ (,s) + f 2 ~ (s) - V{(s)) ds < 0, 

J -72 

anc? /or z € K \ I 2 , we have the analogous inequality: 



(63) / Re (f 2 -(s)+f+( S ) + fi( S )-V 2 '(s)) ds<0, 

J *y 2 

Remark: We note that for the quartic case considered here, Theorem 16 . II solves the Existence Problem 
as described in Section 3, for all a > A. 

Proof of (1) and (2): Straightforward asymptotic calculations for z — > 00 (which are left to the diligent 
reader) show that the three roots {r 3 }| =1 have Laurent expansions of the form (|5~Hl - lt56|l . and it is immediately 
clear that l|58(l holds true. 

Proof of (3) and (4): Whether one calculates the resultant of the polynomial 

P3\z, w) = w — z w + (z + a) w + a z + [3z 

with its derivative (d j dw)ps(z , w) or just solves the equation ij53"|) directly, one arrives at the following 12th 
degree polynomial in the variable z whose roots are potential branch points for the roots {rj(z)}j =1 : 

(64) q(z- a, p, a) = 36a 2 z 12 + 9(4/3 + l)z w + 9 (2a - 18a 2 ) z 8 + 9 (-27a 4 - 18aa 2 + a 2 - 18/3 - 4) z 6 + 

+9 (-54/3a 2 - 12a - 18a/?) z 4 + 9 (-12a 2 - 27^ 2 ) z 2 - 36a 3 . 
The polynomial q is even and calculations are greatly simplified by considering 

(65) q(t) = q{Vt) = 36a 2 t 6 + 9(4/3 + l)t 5 + 9 (2a - 18a 2 ) t A + 9 (-27a 4 - 18aa 2 + a 2 - 18/3 - 4) t 3 + 

+9 (-54/3a 2 - 12a - 18a/3) t 2 + 9 (-12a 2 - 27/3 2 ) t - 36a 3 . 

The goal is to find (a, /3) so that this polynomial has exactly two positive simple roots and two double roots 
(that will turn out to be complex). We consider the resultant of q(t) with its derivative q'(t): 

(66) Resultant [q(t),q'(t)] = B 1 (a,0,a)B 2 (a,f3,a), 
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where 

(67) Bi(a, ft a) = 729aa 6 + 243 (a 2 - 3/3) a 4 + 27a (a 2 - 15/3 - l) a 2 + a 4 + 27/3(3/3 + l) 2 - 36a 2 ft 

(68) B 2 (a, /3, a) = 729aV° + 729 (3a 4 - (3a 2 - /3 3 ) a 8 + 27 (81a 5 + (8 - 171/3)a 3 - 72/3 3 a) a 6 + 

+27 (27a 6 - 20(9/3 + l)a 4 - 2/3 (28/3 2 - 47/3 + 4) a 2 + 8/3 3 (12/3 - 1)) a 4 + 

+ (-27(36/3 + l)a 5 + 4 (216/3 3 + 810/3 2 + 189/3 + 4) a 3 + 576/3 3 (4/3 + l)a) a 2 + 
-/3(4/3 + l) 2 (-27a 4 + 8(9/3 + 2)a 2 + 16/3 2 (4/3 + 1)) . 

Lemma 6.1. There is A so that for all a > A, the following statement is true: For real a and (3, there is an 
isolated zero of B 2 (a, (3), which is a local maximum, occurring at (a*(a), /3*(a)). The quantities a*(a) and 
(3* (a) possess the following expansions: 

(69) a* (a) = a 2 ( -1 + a^ + (±\ a^ 3 + a" 16 / 3 f] a, (a~ 4 / 3 
/3*(a) = a 4 / 3 I 1 - I £ ) a- 4 / 3 + a-^ 3 > 'ft (a" 4 / 3 



a 4/3 ^l-Q)a- 4 / 3 + a- 12 / 3 f> 



Proof 

For the polynomial B 2 (a, /3) to have a local maximum with value 0, it must vanish together with its first 
partial derivatives. To prove this, we first observe that 

i the resultant of B 2 with d a B 2 (eliminating (3) coincides, up to a constant factor, with the resultant 
of B 2 with dpB 2 : 

a 4 

(71) Resultant [B 2 , 3 a B 2 ] (a) = ——Resultant [B 2 , dpB 2 ] (a). 

ii Moreover, there is a unique value of a (called a*) with the asymptotic expansion (|69|l such that 
the above resultants vanish. There are twelve other roots (of varying multiplicities), which remain 
bounded away from this root by a distance of O (a 2 / 3 ) . 

iii In similar fashion, the resultants of B 2 with its first partials, eliminating a instead of (3, are also 
equivalent modulo a constant factor: 

(72) Resultant [B 2 , d a B 2 ] (/?) = -a 4 Resultant [B 2 ,d t3 B 2 ] (/3). ). 

iv Furthermore, there is a unique value of (3 (called /3*) with the asymptotic expansion (|70|l so that these 
two resultants vanish. There are 12 other roots (of varying multiplicities), which remain bounded 
away from this root by a distance of O (a 4 / 3 ) (independent of a). 

The assertions ii and iv above may be verified by computing the large a behavior of these resultants. 
Next one considers the Sylvester matrix associated to Resultant [B 2 , dpB 2 ] (a) (i.e. the matrix whose 
determinant is this resultant) , and shows that for a — a» , there is an eigenvector of the form 

(73) (/3(a0 10 , /3(a*) 9 , •■■,/?("*), 1)' 

with /3(a*) possessing the first few terms of the asymptotic expansion (|70ll . Since /3(a») must be a root of 
the polynomial Resultant [B 2 , dpB 2 ] (/3), it follows that /3(a») = /3*. 

It now follows that B 2 , along with its first partial derivatives, vanishes for (a,/3) = (a*,/?*). 

Again using the fact that a may be taken large, one may verify that for (a,/3) near (a*,/3*), the surface 
defined by B 2 {a,(3) satisfies 



(74) Hess (B 2 ) (a*, ft) = -4782969 a 80/3 (l + O 



a" 4 / 3 

That an A exists so that for all a > A the above holds true is a straightforward exercise in perturbation 
theory. □ 

Remark: From this point on, we shall take a > A, and a — a*, /3 = ft. However, to avoid cumbersome 
notation, we will usually not emphasize this. For example, q(t) will now refer to q(t; a = a*, (3 = ft), etc. 
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A consequence of Lemma 16.11 is that the polynomial q(t) defined in l|65|) possesses two complex double 
roots. One sees this by first carrying out an asymptotic expansion of q(t) for a — > oo, using the asymptotic 
expansions (|69|) and l|7U|) . to see that q(t) possess the following asymptotic expansion: 



(75) 



a' 6 q{a 2 ' 3 () = 36(C - 1) 2 (C 2 + C + I) 2 + 36 a" 4 / 3 (l + ( + C 2 ) (-3 - 6 ( 2 + ( 3 ) + 
-3 a~ 8/3 (-36 + 9 C - 18 C 2 + 51 C 3 - 6 C 4 + C 5 ) + 
+ 18 a~ 12/3 (-2 + C) (1 - C + C 2 ) - 27a~ 16/3 C 



Perturbative calculations then show that for a sufficiently large, q{a 2 ' 3 C) possesses exactly two real roots. 
These are simple and positive (near C = 1). In fact, the reader may verify that the two roots, which we 
define to be j 2 and j 2 (where 71 and 72 are two positive numbers) satisfy: 



(76) 
(77) 



ll = « 2/3 



ll = a 2 ' 3 



-2/3 



M-a 



Therefore q(t) has four complex roots, appearing in complex conjugate pairs. Since q(t) must have at least 
one double root, there must be two double roots. Therefore q(t) factors as follows: 

(78) q(t) = 36a 2 (t- 7l (a) 2 )(t- l2 (a) 2 ) [(t - A.)(* - X)] * , 

with Imag (A*) > 0. 

Now if one solves the cubic l|53|l . one finds 



(79) 

where 

(80) 

(81) 



1 



w 



-R(z) + 3- 5 / 2 v^qJz^) 



1/3 



H(z) 



-fl( Z )+3^ 5 /V-<?( Z 2 ) 
2 



1/3 



R{z) 



H(z) 



27 



-z° + -z J + 



a 2 z 3 



+ (3z 



and the cube roots appearing in (|79|l are all the same choice of branch, but not necessarily the principal 
branch (the three different choices yield different solutions to the cubic (|53(l h Similarly, the square roots 
appearing in (|79|l are all the same choice of branch, not necessarily the principal branch. (However, different 
choices of the branch of the square root do not yield different solutions.) 

Remark: The representation (|79() . although concise, is at present only useful for local analyses since the 
choice of branches, usually a delicate matter, has been avoided. We will use this representation together 
with local analysis near zeros of q(z 2 ) to describe the possible branch points of the solutions to (|53|l . and 
later we will make a particular choice of branches to carry out a global analysis. 



We already know that the only potential branching points where a solution to (|53|l could fail to be analytic 
are the simple roots of q{z 2 ). Yet another asymptotic analysis for large a shows that R(z) does not vanish 
at any of these roots, and so we conclude that there are only four branch points, namely ±71(0) and ±72(0), 
and at each of these points, the solutions to l|53(l exhibit square-root singularities. (So, for example, in a 
vicinity of 71 , the solutions are analytic functions of y/z — 71. ) 

At each of the four branch points, the square root appearing in l|79|) vanishes, while R(z) does not, and 
so a local analysis in a vicinity of each branch point can be obtained by choosing a particular branch for the 
square root (as will be done below), and then, for each distinct choice of branch of the cube-root, use the 



Taylor expansions 

(82) (-i? + 3- 5 / 2 v^(^)) 



1/3 



(83) (-i?-3- 5 / 2 y^^y 



1/3 




where the phase <j> = <j>(z) is the complex argument of —R(z). Different choices of this phase yield the three 
different solutions to (|53|l . 

However, before doing this, it is advantageous to eliminate the radical from the denominator of l|79Jl : 

1/3 



(84) 



H(z) 



-fl( Z )+3- 5 / 2 y^(^) 



1/3 



-2i0/3 



-R(z) - 3- 5 / 2 y/-q{z 2 ) 



Now using the convergent Taylor expansions l|82|) and l|83|) in (|79|l . the reader may verify (with some 
sweat) that exactly one of the three possible branches of the cube-root is analytic in a neighborhood of each 
branch point. Indeed for that one choice of cube- root, all odd powers of the square-root happen to vanish! 
For example, for z real and near 72, it turns out that R(z) > (see Lemma l(j . 21 below 1 ) and hence one must 
choose cf> to be an odd multiple of tt. The three choices —tt, 7r, 3ir will yield all possible solutions to the cubic 
equation J53J). Now using the Taylor expansions and one finds 

1/3 



(85) w{z)-\z 



-R{z) + 3-^ 2 y/^q{^) 



H(z) 



-fl( Z )+3- 5 /y~g[> 2 ) 



1/3 



-R 



E 



2c,- 



,j=o, j even 



V 5/ y~ g (z 2 )' 

R(z) 



cos (0/3) 



E 



2ic., 



R(z) J 



sin (0/3) 



/ j=i, j odd \ 

with Co = 1, Cx = 1/3, and clearly the choice <p = 37r is the only one for which all odd powers of the square 
root vanish, yielding a solution l|79|l to l|53(l which is analytic in a neighborhood of 72, while the remaining 
two are not. 

Remark: Observe that for the other two choices of the angle 4>, namely — tt and 7r, the two corresponding 
roots coincide at z = 72. These two roots are discontinuous across I2 = (71, 72)- It follows that the boundary 
value on one side of I2 of one of these two roots must coincide with the opposing boundary value of the other 
root. Similarly, for the two roots that are not analytic across ii, their boundary values from opposing sides 
coincide, and at the endpoints of the interval Ii, those two roots are equal. 

Global analyticity properties 

In order to describe the global analyticity properties of the solutions to (|53|l . we will make the following 
choice for the branch of the square root appearing in (|79|) : 

(86) ^J-q{z 2 ) := 6ia [(z 2 - K)(z 2 - A*)] yj z - 71(a) \J z - 72(a) ^z + 7i(«) V z + 72(a) 

where now yf ■ is the principal branch. The function appearing in l|86(l is analytic in C\{[— 72, —71] U [71, 72]}- 
Using this choice of branch for \J — q(z 2 ) in (|79|l . along with the fact that R(z) is strictly positive on 
I2 , it follows that the one solution which is analytic in a neighborhood of 72 can be continued to a fixed 
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neighborhood of the interval [71, 72]. This solution is then analytic for z G C\ [—72, — 7i] (of course, it cannot 
be entire, and so it must have a branch cut, which can only be across (—72, — 7i))- This solution shall be 
referred to as r 1 {z; a). 

A similar argument shows that there is another solution, which will be called 7*2(2:; a), that is analytic in 
C\ [71,72]- 

Remark: In order to complete the proof of (3) and (4) of Theorem 16.11 we will show that the solutions 
ri(z;a) and r2(z;a) possess the z — > 00 asymptotic expansions (|56f) and (|55f) . respectively. We will only 
consider the case of ri(z;a), as the proof that the solution T2{z;a) possesses the expansion (|55H is entirely 
similar. 

The starting point for this global analysis is the signature of the functions R{z), H(z), and q(z 2 ) in a 
vicinity of 72(a), and we are once again aided by the fact that a may be taken large. Indeed, we have the 
following lemma. 

Lemma 6.2. There is A > so that for all a > A, the following statements are true. 

• There is 77(a) > 72(a) so that the function R(z) is positive on the interval (0,77), and is negative on 
(ry(a), 00). 

• The function H{z) is strictly positive for all z > 
Proof 

Using the rescaling z 1— » a 1 / 3 ^, we find 

(87) a- 3 iV /3 C)=P(C)+a- 4/3 73(C), 
where 



m «' = ((^J« + y 

Now a* and /3* possess asymptotic expansions (cf. H69|) . l|70ll ') so that the coefficients ^f- and -^fj in l|88(l and 
(|89|l arc uniformly bounded, and 

(90) ^--1. % -1- 

a A art 6 

The reader may verify that the leading order term, p(£), possesses one positive real root which approaches 
3 1 / 3 . Moreover, p(Q is positive to the left of this root, and negative to the right. 

Next we consider £ in a fixed interval of the form [^73 , X] . The function j3(£) is uniformly bounded on this 
interval, and so the function a _3 i?(a 1 / 3 C) can only possess roots in a vicinity of the root of p(C). However, 
this is a simple root of p(C), and the derivative (d/dOa^^Rfa 1 / 3 ^) is strictly negative in a vicinity of this 
root, guaranteeing that a~ 3 R(a 1 / 3 () possesses exactly one simple root in the interval [-^jj, X]. 

Lastly one considers the behavior for £ large, where the function p(() dominates the function p((), and 
one may then choose X large enough (but fixed) so that ^(a 1 / 3 ^) possesses no roots beyond X, and this 
finishes the proof of the first claim of the Lemma. 

The second claim can be established by similar reasoning, using the asymptotic representation 

a^Hia^O^Uc 6 -*^))-^ 3 (U 2 



9 V \a 2 JJ \3 
It is straightforward to verify that the leading order term on the right hand side is strictly positive for 
C > 72/a 1 ^ 3 j an d the remaining details of the proof of the second claim are left to the reader. □. 



Remark: For the sake of completeness, we note that the function \J — q(z 2 ) is real and negative for 
z G (71(a), 72(a)), and positive imaginary for z > 72. 

Using Lemma 16.21 it follows that as z traverses the real axis, from 72 to +00, the quantity —R(z) + 
3~ 5 / 2 -\/— q(z 2 ) starts out on the negative real axis, enters the second quadrant, exiting the second quadrant 
by crossing the positive imaginary axis, and then diverges to 00 in the first quadrant. 
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In the text following (JSSJ, we have shown that for z € M and near 72, the appropriate choice of cube root 
yielding a function that is analytic is determined by (f> = 3n. As z traverses the real axis from 72 to +00, this 
angle must converge to 2ir, and hence we have determined the behavior of the root r\{£) to be as follows: 



(91) r 1 (z) = iz 3 + i 2 3 e 2 -/ 3 1 + ^ + 



\/3 3 m ( 1 



2z 



O 



4 1 ~ \ ^5 



2 

Z \ Z z 



which establishes <|56ll . 

We will leave to the motivated reader the proof that the root r 2 , which is analytic in C \ [71 , 72] , possesses 
the asymptotic expansion (|55|l as z — ► 00. 

This completes the proof of (3) and (4) of Theorem 16. II 

Proof of (5) of Theorem 16.11 To see that the boundary relation (|59|l holds true, one first observes 
that for z G I\, both boundary values rf are roots of the cubic (|53|) . Since they are not both the same root, 
and since r 2 is analytic in a neighborhood of I\ , it follows that there must be a relation between rf and rf . 
The reader may verify that this monodromy relation is 



f for z £ I\. 



It is then straightforward to verify that 

(fi + fi) + h=rt + rf + r 2 +a = rf + r 2 + r+ + a = z 3 + a, 

since the sum of the roots of the cubic l|53|l must be z 3 . The boundary relation H6U|) follows by very similar 
reasoning, the difference being that on I 2 , the boundary values rf are related to the boundary values of r% 
as follows: rf = . 

Proof of (6): The inequalities l|61f) may be seen to be true by one of the following two approaches. We 
will consider j = 1, as the considerations for j = 2 are entirely similar. 

(1) Using the expansions (|82|l and (|83|) (with suitable choice of the angle <fi), valid over the entire interval 
II, one finds that the difference between the boundary values from the upper and lower half planes 
is positive imaginary. 

(2) One shows that the two roots rf and rf are complex conjugates of each other, and that rf has 
positive imaginary part. 

Proof of (7): We now prove that the inequality (|62|l is true. First, recall that on (|59|l holds true. 
Next observe that for s £ R \ Ji, 

(92) /f( S ) + f+ {a) + f 2 -(s) - V{(s) = 2h{s) + / 2 "( S ) - s 3 - a = n(s) - r^(s). 

It turns out that one must verify the inequality (|62() separately on each of the intervals (—00, —72), (—71, 71), 
h, and (72,00). 

Local analysis near the branch points —72, —71, and 72 show that on each of the intervals (—00,-72), 
(—71,71) and (72,00) there are three real roots. Since the only places where any two roots can coincide 
are the branch points ±71 and ±72, the quantity (|92|l is of one sign for s in each of (—00, —72), (— 71, 71) 
and (72, 00). The signature in each of these intervals can be determined by a (by now) straightforward local 
analysis in a vicinity of each branch point. 

For example, for s = —71, the two roots r\ and rz coincide, and the following local expansions hold true: 
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(93) 



(94) 



-R 



;(2n/3)(l + 0(\q(s 2 )\)) 



2z p-^y^F) 
3 I i?(s) 



3 in(27r/3)(l + 0(| g ( S 2 )|)) 



r 3(s) - -S 3 



-1? 



cos(4./3)(l + 0(| 9 ( S 2 )|)) - | ( 3 " 5/ %^ 



(47r/3)(l + C(|g( 5 2 )|)) 



Now for s real, s > —71, and s near —71, y— q(s 2 ) is a number on the negative imaginary axis. Therefore 

~4i ( z- 5 ' 2 ^qjf) y 



(95) 



n(s) - r 3 (s) 



-R 



R(s) 



sin(27r/3) + C(|g(s 2 )|) < 0, 



and then it follows that the inequality (|62[) is true for s 6 (—71,71)- Similar analysis shows that the same 
inequality is true on (—00, —72), and we will leave the details of those calculations to the reader. 

For z in the interval I2, the root r\ remains real, but the root r 3 (either the "+" or the "— " boundary 

value) is not real, and the quantity Re \r\{z) — \z)j may change signs without two roots coinciding. 

However, recalling the form (|79[) of the roots of the cubic l|53|) . along with the useful relationship l|84|) . we 
may represent the roots 77 as follows: 



(96) 



r, = rz + 



-R(z) + 3- 5 / 2 ^q(z2) 



1/3 



-2i<f> 



,. /3 / -fl(^-3-^V=gg) 



1/3 



where the choice of cube- root, and the phase <f>j, depends upon which root is being considered. 
We will need the following properties: 

• The function yj — q(z 2 ) is real and negative for z S (71(a), 72(a)). 

• The function R{z) is strictly positive for z £ (71(a), 72(a)). 

• The quantities — R(z) + 3~ 5 ^ 2 y // — q(z 2 ) and — R(z) — 3~ 5//2 -\/— q(z 2 ) are nonzero for z € [71(a), 72(a)], 
and consequently they are both strictly negative for all z g [71(a), 72(a)]. 

With the above properties, and the choice 4>\ = 37r (see the text following JSSJ), we have 



(97) 



[n - r 3 



,103/3 



e »T _ e »03/3 



-R(z) 4 


- 3" 5 / 2 y/-q{sP) 




2 


-R(z)^ 


-i- 5/2 ^J-q{z 2 ) 




2 


-R(z) 


+■ 3- 5 / V-g(z 2 ) 


2 



1/3 



-i?(z)-3- 5 / 2 v^M 



1/3 



1/3 



+ e 



1/3 



-2i0i/3 e i0 3 73 



-103/3 



-i?(z)-3- 5 / 2 V-g(^ 2 ) 



X/3N 



- 3" 5 /V-?( z2 ) 



1/3 



and therefore we have shown that Re (r\ — r^) < for all z G [71, 72], which completes the proof that (|62l) 
is true for z € 71 , 72 ) • 

Lastly, since Re m — r^) < for 2 € [71,72], an d since all three roots are distinct for all z € (72,00), 
it follows that this inequality remains true in (72,00) as well, and hence we have completed the proof that 
(JSU holds true for all 2 6l\/i. 

The proof that l|63|) holds true follows by very similar arguments, and we will omit these details here. 

Having completed the proof of (7) of Theorem 16. II we have completed the proof of the entire Theorem. 

□ 
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Remark: We note that at each of the four endpoints, the roots r\ and r 2 possess square root singularities. 
This implies that in addition to verifying all of the equalities and inequalities of the ideal situation A-D, 
conditions El and E2 also hold true. 



7. Complete analysis of the Riemann-Hilbert problem II. II in the "one-cut" case 

The starting point of this section is the Riemann-Hilbert problem (|12fl . In this section we will show that 
if the ideal situation described in sectionals achieved (including the "Airy conditions" El and E2), and if 
each of the sets Ij ,j = 1, 2 is a single interval, then one may obtain a complete asymptotic description of the 
solution of the Riemann-Hilbert problem l|12fl . Blcher and Kuijlaars [3] have carried out a Riemann-Hilbert 
analysis for the Gaussian case, and their work may be applied directly to the present more general situation. 
Therefore, we will only present a summary of the relevant transformations. Of course, it should be possible 
to generalize the work of this section to the case of arbitrary (finite) numbers of subintervals comprising the 
sets I\ and I 2 , but this is putting the cart before the horse as there is at present no existence theorem in the 
multi-cut case. In what follows we will let the endpoints of the intervals be defined as follows: 

(98) 7 1 = (7iA), h = {12,82), 

and we will use the following particular choices of antiderivatives of /1 and f 2 to define g\ and <7 2 : 

(99) .91 = — / h(s)ds- Cl , 

ni J Sl 

(100) 92 = — f f2(s)ds-c 2 , 

n-2 Js 2 

where in (|99|) the contour is chosen in the cut plane C \ (— 00, Si), and in (|100f> the contour is taken in 
C \ (—00, 62), and the constants c\ and c 2 are chosen so that © holds true, namely 

(101) Cl = ilog5 1+ ^ + ^( S )-1^ S) 

(102) C2 = \\0g82 (f 2 (s) + ^) ds. 

(Of course, then go is determined via (|10|) . % ) 

Remark: In what follows, we will assume a certain type of genericity in the following sense. We are 
considering a sequence of triples, (m, n 2 , n) so that n = ?ii+?i 2 and rij/n — > Xj. The implicit characterization 
of the functions gi and <? 2 in terms of the equation E(w, z) = with E(w, z) defined in 1)36(1 depends on x\ 
and X2, and we will assume that the "ideal situation" holds true not just for x\ and £ 2 , but for all values of 
{x'^x'2) near (ii,a; 2 ). 

The starting point is to realize that, because all of the conditions of the "ideal situation" are satisfied, 
the jump matrix appearing in (|12|l . namely 

( e - n (9o -So) e ng ~+n 1 g+-nV 1 -n 1 e 1 ^ng' +n 2 g} -nV 2 -n 2 l2 
e nx{gt-9T) 

\ e n 2(g}-g 2 ) 



(103) V B 



takes on one of four forms, as follows. 

(I e —nPi{x) e -nP 2 { 

1 I' iorxe (-00,71) U(<5i, 72 )U(<5 2 , 00), 

1 J 

( e -in 2 8 2 (x) g-nPi(i) ^ \ 

1 , for x e h 

e ^ 2 e 2 (x)J 

/g-miflifi) ^ e -nP 2 (x)\ 

(106) V B (x) = einiOiix) , for a; eli, 



where 



• Pi = Xx(gf + g 1 ) + 2^2.92 — Vi — x\i\ is strictly negative on (—00,71) U (Si, 00), 

• Pi = x 2 {g2 + 9i) + ^151 —V2 — £2^2 is strictly negative on (-00, 72) U (5 2 , 00), 

• #2(2:) = —iigt ( x ) ~92~( x )) i s strictly positive and strictly decreasing on I2, and possesses an analytic 
continuation (also called 82) off the interval I2, 

• 9i(x) = —i(g^(x) — gi(x)) is strictly positive and strictly decreasing on Ii, and possesses an analytic 
continuation (also called 9i) off the interval Ii, 

7.1. Second transformation: Opening lenses. Figure 1 below shows a decomposition of the plane into 
five regions, formed by defining four lens shaped regions surrounding the intervals Ii and 12- We will define 
the matrix D in each of these regions, as shown in Figure 1. (In the region exterior to the four lens shaped 
regions, D = B.) The matrix D is clearly piecewise analytic. Because it is related directly to the matrix B, 
it also solves a Riemann-Hilbert problem. However, we will not write this problem down, as it is somewhat 
auxiliary to our final goal. This definition is quite similar to the definitions (5.2), (5.4) and (5.6) of [3J- 





0> 


Figure 1. The decomposition of the plane into 5 regions, the interior of the 4 lens shaped regions 
surrounding the intervals (71, Si) and (72, S2), and the one exterior region. The matrix D is defined 
as shown in each of the 4 bounded regions, and D = B for z exterior to all four lens shaped regions. 

Remark: In order to define the matrix D, one must extend, for j = 1,2, the quantities — gj and 
Pj, for j = 1,2 off the interval The justification for this uses the connection to the algebraic curve. For 
example, 



(107) 



1 

X2 



/+ - f 2 ~dx 



fi 



fS 



dx 



fo 



/a 



dx, 



which clearly demonstrates the local analyticity of g£ — g 2 ■ The other quantities are seen to be analytic by 
similar calculations. 

The next step is to define a piecewise analytic matrix valued function D approx (z) which, it will be shown, 
is a globally uniform approximation to the matrix valued function D. We will require the following auxiliary 
functions. Set 



(108) 
(109) 
(110) 



P 



P, 



^2(72), 



Ps 2 



n(Si), 

r2(S 2 )- 



Let r+ denote the image of the +-side of the interval Ii under the transformation £(z) = ro(z) (recall that 
tq(z) has branch cuts along Ii and ^2), and let T2 denote the image of the +- side of the interval I2 under the 
same transformation £(z) = ro(z). The reader may verify that for each j = 1,2, is an arc in the upper 
half plane connecting P lj and Psj . 
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Now we define three functions, Mj{z) (j = 1,2,3) as follows (cf. (6.9) of [3])- 

e-a 2 



(111) M x (0 = 



(112) M 2 (0 = c 2 

(113) M 3 (0=c 3 



y/G - P-n ) (£ - Ps 1 ) (£ - ^» ) (I - ^ 2 ) ' 



V(C-^i)(e--P«i)(e--P7 2 )K-ft 2 ) : 



where in each instance, the quantity — P 7l ) (£ — -P^ ) (£ — -P 72 ) (£ — -Pfo ) i s taken to be analytic in C \ 
(rj" U rjj , and the constants ci and c 2 are taken to be 



(114) a = ^y/(a-P m )(a- P Sl ) (a - P 72 ) (a - Pg 2 ), 



(115) c 2 = ^V(-a - iV) (-a - ^) (~ a ~ P i*) (~ a ~ P ^)- 
We will require the following 3x3 matrix P(z): 

/Mi(r (z)) A/iMz)) Mx(r 2 (z)) N 

(116) P(z) = M 3 (r (z)) M 3 (n(z)) M 3 (r 2 (z)) 

\M 2 {r (z)) M 2 (n(z)) M 2 {r 2 {z)) / 

For future reference, the matrix P(z) solves the following Ricmann-Hilbcrt problem. 

Riemann-Hilbert Problem 7.1. (a) P is analytic onC\ (Ii U I 2 ). 
(b) The boundary values of P satisfy 

(117) P+(x)=P.(o 






1 
































1 














(118) P+(x)=P-(i 
(119) 

(c) As z — > 00, we have 
(120) P(z)=(l + o(~ 

Next we draw 4 small circles, centered at each of the endpoints 71,61,72,62, and decompose the plane 
again, as shown in Figure 2. 







7i / V 01 J \J2 



Figure 2. Four disks, centered at each of the endpoints of the intervals l\ and I 2 . We decompose 
the plane into five regions: the four disks and the one exterior region. 

We will now define a matrix D approx which will be a uniformly valid approximation to the matrix D. The 
matrix D approx will be defined separately in each of the 5 regions as shown in Figure 2. We set 

(121) D approx (z) = P(z), for z outside the four disks. 

For z in the interior of each of the four disks, D approx is defined in a (by now) standard way using the 
solutions of the Airy equation. Rather than presenting the details for each of the disks, we will focus only 
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on the disk centered at 8± , as the construction for each of the other disks may be carried out by very similar 
calculations. Define 



(122) 



S(z) 



61 



(2f 1 (x)+f 2 (x)-V;(x))dx 



-1 2/3 



taken to be analytic in a neighborhood of z = Si, and with S'(S\) > (which is possible because of 
conditions El and E2). It turns out that with this definition, one may take an appropriate branch so that 
S(z) 3 / 2 = — x^i( z ) f° r z m t nc upper lens shaped region above I\, and S(z) 3 / 2 — ^-6i(z) for z in the lower 
lens shaped region below I\ . 

Next, define the matrix as follows (cf. gj (7.13)]). 

o 1, 

-y'i(s) -y 2 (s) 
01, 

'-y 2 {s) yi (s) 0\ 

1, 

^o(s) 2/1 («) ON 

2/o(«) 2/i («) 
1, 



(123) 



9(s) 



for < args < ?f, 



for ^ < arg s < n, 



for — 7r < args < — 



for 



< arg s < 0. 



where ?/0j ?/ij and ?/2 are the following solutions of the Airy equation (recall that Ai(x) is the unique solution 
to Airy's equation y" = xy which is exponentially decaying for x — > +00): 



(124) 



y (s) = Ai(s), j/i(s) = wAi(ws), y 2 (s) = w 2 Ai(w 2 s) 



(here w = e 3 ). 

And then define E n (z) for 2 in the disk centered at 5i (again following as follows: 



(125) 



1 -1 
E n (z) = yfrP(z) ( -i -i 
1 




0^ 

n-V*S(z)-V* 
1, 



Finally then, we may define D approx (z) for z in the disk centered at 5\\ 



(126) 



(z) = E n (z)^>(n 2 ^S(z)^ diag(e-^(^-^), e^H 9 ^ ), l). 



' approx 

There is a similar definition of D a p prox in each of the other 3 disks shown in Figure 2, which are entirely 
similar to the construction presented here j3j- We will omit these definitions here. 



The matrix D 



approx 



is a piecewise analytic function in the entire plane, and it is a global approximation 
to the matrix D(z) defined in Figure 1. The proof of this follows a by now straightforward procedure, which 
has been described in (8.1)-(8.6)]. Indeed, setting 



(127) 



R(z) = D(z)D- 1 prox (z) 



one verifies that R(z) satisfies a Riemann-Hilbert problem of exactly the same type as the analogous function 
R(z) defined in [3J (8.1)]: jumps across the boundaries of the disks are I + C(n _1 ), and jumps across all 
other contours are I + O (e~ cn ), and in addition the actual contours chosen may be deformed slightly (as in 
the text between (8.5) and (8.6) of |3j). Therefore one may conclude that 



(128) 

uniformly for z£C. 



R{z) = 1 + 



n(\z\ 



00, 
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8. ASYMPTOTICS FOR EIGENVALUE STATISTICS OF RANDOM MATRICES WITH SOURCE 



In this short section we will describe some of the results that follow immediately from the results of the 
previous section. The proofs of these results follow by considerations entirely similar to those carried out 
in the Gaussian case by Bleher and Kuijlaars [3]: one expresses the kernel K n (x,y) in terms of the explicit 
transformations leading to R(z), whose asymptotic expansion for n — > oo is under control. 

Theorem 8.1. Suppose that the conditions described as the "ideal situation" ' , items A-D, described in Section 
2 are satisfied, along with the conditions described in El and E2 of Section 3. Suppose in addition that the 
sets I\ and I2 are each single intervals, with I\ = (71, £1) and I2 = (72,82)- Then the following results hold 
true. 

• The mean density of states, p n {x) := K n (x,x) converges, as n — > oo, to a limiting mean density, 
supported on I\ U I2, with 

(129) Urn Pn (x) = (sH*) - 97 W) > f or x ^ ^ 3 = h 2. 

n— >oo 2lTl dx J J 

• Bulk universality holds true: for every xq € (71, (Si) U (72,^2)7 omd every we have 

1 / u v \ shnr(u — v) 

(130) hm — - — -K n x H 7 — r,x Q H — - = r— . 

n-too np(xo) \ np(xo) np(xo) / ir(u — v) 

• Edge universality also holds true, at each endpoint of I\ U I2 •' There are constants c 7l , cs ± , c 72 , and 
cg 2 so that for j = 1, 2, the following statements hold true: for every we have 

I f u v \ _ Ai(u)Ai'(v) - Ai(v)Ai'(u) 

( 131 > ^ ( C7jrl )2/3 A « [V + Jc~n^ nj + (c 7j n)2/3 J - —v ' 

1 / u v \ Ailu)Ai'(v) - Ai(v)Ai'(u) 

( 132 ) 7—3273^ S i 7— ^273 . S > 



( % n) 2 /3 V 3 (%n)2/3' J (c S] n)V3 / 

Remark: As proven in Sectional all of the conditions A-D along with El and E2 hold true for all a 
sufficiently large, for the quartic case V(x) = jx 4 , with n even and ?ii = 712 = n/2. Thus bulk and edge 
universality have been extended beyond the Gaussian case for random matrices with source. 

It is natural to expect that the same results should hold true for all a sufficiently large, under the more 
general assumption that the external field V(x) is convex, real analytic, and with sufficient growth for 
\x\ — > 00. While it seems a rather daunting task to carry out the requisite asymptotic analysis for a —> 00 to 
prove that a suitable curve exists (note that it need not be algebraic), there is a possibility to combine the 
analysis described herein with an analysis of the coupled variational problem described following formula 
(1.15) in 0. 
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